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ABSTRACT 



To model the dynamics of the chiral order parameter in a far from equilibrium 
phase transition, we consider quenching in the 0(4) linear sigma model. We ar- 
gue, and present numerical evidence, that in the period immediately following the 
quench long wavelength modes of the pion field are amplified. This results in large 
regions of coherent pion oscillations, and could lead to dramatic phenomenological 
consequences in heavy ion collisions. 
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Among the most interesting speculations regarding ultra-high energy hadronic 
or heavy nucleus collisions is the idea that regions of misaligned vacuum might 
occur [1,2,3,4,5]. In such misaligned regions, which are analogous to misaligned 
domains in a ferromagnet, the chiral condensate points in a different direction 
from that favored in the ground state. If we parametrize the condensate using the 
usual variables of the sigma model, misaligned vacuum regions are places where 
the four-component (a, 7?) field, that in the ground state takes the value (v, 0), is 
partially aligned in the n directions. If they were produced, misaligned vacuum 
regions plausibly would behave as "pion lasers", relaxing to the ground state by 
coherent pion emission. They would produce clusters of pions bunched in rapidity 
with highly non-Gaussian charge distributions. Thus a misaligned vacuum region 
starting with the field in the n\ — 112 plane would emit only charged pions (equally 
positive and negative, since the fields are real), while a misaligned vacuum region 
starting with the field pointing in the 1T3 direction would emit only neutral pions. 
More generally if we define 

R= ^ (1) 

then assuming that all initial values on the 3-sphere are equally probable, the 
probability distribution V(R) is given by 

V(R) = l -R- 1 ' 2 . (2) 

As an example of (2), we note that the probability that the neutral pion frac- 
tion R is less than .01 is 0.1! This is a graphic illustration of how different (2) is 
from what one would expect if individual pions were independently emitted with 
no isospin correlations. According to (2) there is a substantial probability that, 
say, a cluster of 70 pions would all be charged - something that could essentially 
never occur for incoherent emission. There may be some hint of such behavior in 
the Centauro (overwhelmingly charged) and anti-Centauro (overwhelmingly neu- 
tral) events reported in the cosmic ray literature. [6,7] Here we propose a concrete 
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mechanism by which such phenomena might arise in heavy ion collisions, for which 
the plasma is far from thermal equilibrium. 

In previous work [5,8] we considered the equilibrium phase structure of QCD. 
We argued that QCD with two massless quark flavors probably undergoes a second- 
order chiral phase transition. For many purposes it is a good approximation to 
treat the u and d quarks as approximately massless, and we expect that real QCD 
has a smooth but perhaps rapid transition as a function of decreasing temper- 
ature from small intrinsic to large spontaneous chiral symmetry breaking. At 
first sight it might appear that a second-order phase transition is especially fa- 
vorable for the development of large regions of misaligned vacuum. Indeed the 
long-lived, long- wavelength critical fluctuations which provide the classic signature 
of a second-order transition are such regions. Unfortunately the effect of the light 
quark masses, even though they are formally much smaller than intrinsic QCD 
scales, spoil this possibility [5]. The pion masses, or more precisely the inverse 
correlation length in the pion channel, are almost certainly not small compared to 
the transition temperature. As a result the misaligned regions are modest affairs at 
best even near the critical temperature. They extend at most over a few thermal 
wavelengths, and almost certainly do not contain sufficient energy to radiate large 
numbers of correlated pions. 

Here we will consider an idealization that is in some ways opposite to that of 
thermal equilibrium, that is the occurence of a sudden quench from high to low 
temperatures, in which the (a, 7?) fields are suddenly removed from contact with a 
high temperature heat bath and subsequently evolve mechanically. We shall show 
that long wavelength fluctuations of the light fields (the pions in QCD, which would 
be massless if not for the quark masses) can develop following a quench from some 
temperature T > T c to T = 0. The long wavelength modes of the pion fields 
are unstable and grow relative to the short wavelength modes. 

Before we enter into the details of our model and simulations it seems appropri- 
ate to discuss the qualitative reason for the result just mentioned, which we suspect 
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may be of wider interest. Whenever one has spontaneous breaking of a continuous 
global symmetry, massless Nambu-Goldstone bosons occur. The masslessness of 
these modes occurs through a cancellation following the schema 

m 2 = -fi 2 + A0 2 , (3) 

where the second term arises from interaction with a condensate whose expectation 
value ((p) 2 satisfies ((f)) 2 = (J?/\ = v 2 in the ground state. However following a 
quench the condensate starts with its average at rather than v, and generally os- 
cillates before settling to its final value. Whenever ((f)) 2 < v 2 , m 2 will be negative, 
and sufficiently long wavelengths fluctuations in the Nambu-Goldstone boson field 
will grow exponentially. The same mechanism will work, though less efficiently, if 
one has only an approximate symmetry and approximate Nambu-Goldstone modes, 
or for that matter other modes whose effective m 2 is 'accidentally' pumped negative 
by their interactions with the condensate. 

Returning to QCD, we shall use the classic Gell-Mann-Levy Lagrangian [9] to 
describe the low energy interactions of pions: 

C= J d 4 x{ X - d^ a d i( f) a - ^(« Q - v 2 ) 2 + Ha) , (4) 

where is a four-component vector in internal space with entries (a, tt). Here, A, 
v, and H oc m q are to be thought of as parameters in the low energy effective 
theory obtained after integrating out heavy degrees of freedom. We shall treat (4) 
as it stands as a classical field theory, since the phenomenon we are attempting 
to address is basically classical and because as a practical matter it would be 
prohibitively difficult to do better. We shall be dealing with temperatures of order 
200 MeV or less, so that neglect of heavier fields seems reasonable. 

To model a quench, we begin at a temperature well above T c . The typical con- 
figurations have short correlation lengths and ((f)) ~ 0. (((f)) ^ because H ^ 0.) 
One then takes the temperature instantaneously to zero. The equilibrium config- 
uration is an ordered state with the (f> field aligned in the a direction throughout 
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space, but this is not the configuration in which the system finds itself. The actual, 
disordered configuration then evolves according to the zero temperature equations 
of motion obtained by varying (4). Quenching in magnet models has been much 
studied in condensed matter physics [10]. However in that context it is usually 
appropriate to use diffusive equations of motion, because the magnet is always 
in significant contact with other light modes (e.g. phonons). For this reason the 
condensed matter literature we are aware of does not directly apply to our problem. 

Let us now consider how the physics of quenching may be applied in relativis- 
tic heavy ion collisions. In Bjorken's [11] picture of such a collision, the incident 
nuclei as seen in the center of mass frame are both Lorentz contracted into pancake 
shapes. They pass through each other, and leave behind a region of hot vacuum. 
The baryon number of the incident nuclei ends up in that part of the plasma head- 
ing approximately down the beam pipe, and the central rapidity region which we 
consider in this paper consists of plasma with approximately zero baryon num- 
ber which expands and cools through T = T c and eventually hadronizes into the 
detected pions. If the collision is energetic enough to create a region of plasma 
well above T = T c , the field will indeed be fluctuating among an ensemble of 
disordered configurations. The plasma cools rapidly as it expands, and cannot be 
exactly in thermal equilibrium. If it cools fast enough, the configuration of the <fi 
field will "lag," and as in a quench the system will find itself more disordered than 
the equilibrium configuration appropriate to the current temperature. 

We have done numerical simulations of quenching to zero temperature in the 
linear sigma model with an explicit symmetry breaking term Ha which makes 
the pions massive. Turok and Spergel [12] have considered this scenario with no 
explicit symmetry breaking term as a cosmological model for large scale structure 
formation in the early universe. They find a scaling solution in which the size of 
correlated domains grows at the speed of light as larger and larger regions come 
into causal contact and align. When the O(N) symmetry is explicitly broken, 
however, we do not expect a scaling solution. The Ha term tilts the potential, the 
vacuum manifold is not degenerate, and in a time of order m" 1 the scalar field in 
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all regions (whether in causal contact or not) will be oscillating about the sigma 
direction. 

In numerically simulating a quench, we choose and randomly independently 
on each site of a cubic lattice. This means that the lattice spacing a represents 
the correlation length in the initial conditions. If the initial conditions are cho- 
sen to model a thermal ensemble at some temperature T > T c , then the O(N) 
symmetry is not spontaneously broken, and the lattice spacing a represents the 7r 
and a correlation lengths which are approximately degenerate. Because the initial 
conditions must not distinguish between the a and n directions in internal space, 
we must use the linear sigma model (4) instead of integrating out the a which is 
heavy at zero temperature to obtain the nonlinear sigma model. In the simulation 
whose results are shown in Figure 1, we chose and randomly from gaussian 
distributions centered around = = and with root mean square variance v/2 
and v respectively. The three parameters v, H, and A in (4) determine m n , m a , 
and f n = (0|er|0) according to 

A(0H0)((0H0) 2 - v 2 ^j - H = , (5) 

m " = WW) ' and m ° = 3A <°H°> 2 -^ 2 • ( 6 ) 

Note that f-^ = (0|<r|0) > v for H ^ 0. In interpreting our results we must 
remember that while in the code we are free to choose the energy scale by setting 
the lattice spacing a = 1, a actually represents the initial correlation length. In 
choosing the parameters for Figure 1 we assumed that a = (200 MeV) -1 and then 
chose v = 87.4 MeV = 0.4372 a -1 , H = (119 MeV) 3 = 0.2107 a~ 3 , and A = 20.0 
so that f n = 92.5 MeV, m n = 135 MeV, and m a = 600 MeV. 

With parameters and initial conditions chosen, we evolve the initial config- 
uration according to the equations of motion using a standard finite difference, 
staggered leapfrog scheme. (For details of the numerical method see [13,14] and in 
particular equations 32-35 of [14].) We used a time step dt = min(^, jjf^-)- We 
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verified that our results do not change if the time step is reduced. After each two 
time steps, we computed the spatial fourier transform of the configuration, and 
from that obtained the angular averaged power spectrum. In Figure 1, we plot the 
power in modes of the pion and sigma fields with spatial wave vectors of several 
different magnitudes k as a function of time. All the curves start at approximately 
the same value at t — because the initial power spectrum is white since we chose 
independently at each lattice site. The behaviour of the low momentum pion 
modes is striking. While the initial power spectrum is white and, as ergodicity 
arguments would predict, the system at late times is approaching an equilibrium 
configuration in which the equipartition theorem holds, at intermediate times of 
order several times m~ x the low momentum pion modes are oscillating coherently 
in phase with large amplitudes. 

The simulation was done in a 64 3 box. We verified that finite size effects are 
not important even for the longest wavelength (ka = 0.20) mode shown in Figure 1 
by checking that the behaviour of the ka = 0.31 mode is the same in a 32 3 and 
a 64 3 box. The behaviour of only one component of the pion field is shown in 
Figure la. The other two look qualitatively the same. It should be noted that the 
exact height of the peaks in the curves depend on the specific initial conditions. 
If the simulation is run with the same initial distribution for <fi and <fi but with a 
different seed for the random number generator, the heights of the peaks change, 
and the relative sizes of the peaks in the three different pion directions change. 
The qualitative features of Figure 1 — the growth of long wavelength modes of the 
pion field — do not depend on the specific realization of the initial conditions. We 
discuss below how the results change as a function of parameters in the potential 
and the initial distributions. 

Let us compare the growth of long wavelength modes found numerically with 
expectations from the mechanism previously discussed. Suppose the potential V 
in (4) were simply V(<f>) = (m 2 /2)(f> a (f> a . Then, the equations of motion would 
be linear, and modes with different spatial wave vector k would be uncoupled. 
Each curve in Figure 1 would be a sinusoid with period tx j \/ m 2 + k 2 and constant 
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amplitude. (The power spectrum, being quadratic in the fields, oscillates with one 
half the period of the fields.) The period of the oscillations in Figure la is indeed 
given by -k / \j m 2 + k 2 , but the amplitudes are far from constant. This behaviour 
can be qualitatively understood by approximating 2 in the nonlinear term in the 
equation of motion by its spatial average: 

Q Q (f,t)~(0 Q Q )(t) . (7) 

This approximation is exact in the large N limit [12]. In our problem, N = 4. 
Using (7) and doing the spatial fourier transform, the equation of motion for the 
pion field becomes 

d 2 

-^7f(k,t) = -m 2 eff (k,t)n(k,t) (8) 
with the time dependent "mass" given by 

m 2 eff (k, t) = -\v 2 + k 2 + A(0 2 )(t) (9) 

where k = \k\. 

Figure 2 shows the time evolution of (0 2 ) in the same simulation whose results 
are shown in Figure 1. In the initial conditions, is gaussian distributed with 
(<p 2 ) < v 2 . Therefore, for a range of wavevectors with k less than some critical 
value, rn 2 ff < and the long wavelength modes of the pion field start growing 
exponentially. (0 2 ) grows and then executes damped oscillations about its ground 
state value (0|<r|0) 2 . A wave vector k mode of the pion field is unstable and grows 
exponentially whenever (0 2 ) < v 2 — k 2 /\. Modes with k 2 > \v 2 can never be 
unstable. The k = mode is unstable during the periods of time when the (0 2 ) 
curve in Figure 2 is below v 2 . Since (0 2 ) is oscillating about (0|<r|0) 2 > v 2 , after 
some time the oscillations have damped enough that (0 2 ) never drops below v 2 , 
and from that time on all modes are always stable and oscillatory. In general, 
longer wavelength modes are unstable for more and for longer intervals of time 
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than shorter wavelength modes as (<fi ) oscillates. Also, "^f// * s more negative 
and more growth occurs for modes with smaller k. Making the approximation (7) 
and thus using (8) cannot be expected to completely reproduce the effects of the 
nonlinear term in (4) which is local in position space. Nevertheless, it predicts 
that the long wavelength modes of the pion field go through alternating periods 
of oscillatory behaviour and exponential growth spurts and therefore gives us a 
good understanding of the behaviour of these modes in Figure la. Because the 
timing of the growth spurts for different (low) k modes are all determined by 
(0 2 ), different modes have their growth spurts at the same times and all the long 
wavelength modes in Figure la oscillate in phase. The longer wavelength modes 
have more, stronger, and longer growth spurts, and are therefore amplified more 
as in Figure la. In an equilibrium phase transition, explicit symmetry breaking 
keeps the correlation length at T c finite and, in the case of QCD, too short to be of 
interest. [5] After a quench, on the other hand, arbitrarily long wavelength modes 
of the pion field are amplified even though the pion mass is non-zero. 

What about the a modes? For a 2 < t> 2 /3, the effective mass 2 for the k = 
mode is negative. This condition is far less likely to be satisfied than (0 2 ) < v 2 . 
Hence, while some growth of the low momentum a modes is possible, as Figure 1 
confirms, they do not grow as much as the low momentum pion modes do. The 
behaviour of the sigma modes is also complicated by the fact that oscillations with 
two periods — 27i/m n and 2n/m a — contribute. Roughly, this occurs because 
oscillations ostensibly along the pion directions are oscillations in a curved valley 
and are therefore seen in sigma also. 

At late times ((f) 2 ) oscillates with small enough amplitude that it never goes 
below v 2 and consequently no modes are ever unstable. Hence, if we make the ap- 
proximation (7) we would expect that at late times each mode in Figure la would 
continue to oscillate with approximately constant amplitude, with the longer wave- 
length modes maintaining the large amplitudes acquired during their exponential 
growth spurts. That this is not what is seen in Figure la reflects the effects we 
neglected in making the approximation (7). Because the modes are in fact cou- 
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pled and the equations of motion are actually nonlinear, ergodicity arguments 
suggest that eventually equipartition should apply. If the energy is equally divided 
among modes, then the (amplitude) 2 in a mode should be inversely proportional 
to (m 2 + k 2 ). Also, the amplitude of the sigma modes at late times should be less 
than that of pion modes of the same k because m a > m n . The results shown in 
Figure 1 are consistent with the assumption that at late times equipartition ap- 
plies, although at t — 80a the longest wavelength pion modes are still decreasing in 
amplitude. (We have verified that at even later times they do reach equipartition.) 
It is reasonable that longer wavelength modes take longer to decrease in ampli- 
tude than shorter wavelength modes. As the system heads toward equipartition, 
power has to flow from low momentum modes to higher momentum modes. This 
plausibly happens more quickly at higher momentum because higher momentum 
modes are coupled to a larger number of modes, or, in the continuum, a larger 
volume of momentum space. Hence, we end up with the striking phenomena of 

— # 

Figure la: while the system begins and ends with small Tr(k,t), in between the 
long wavelength pion modes are dramatically amplified. 

In modelling a heavy ion collision, we should include the effects of the expansion 
of the plasma after a quench. In an equilibrium phase transition, the expansion 
causes a decrease in the temperature. After a quench, however, T = and the 
expansion has the effect of reducing the energy per unit volume in the system. The 
simplest way to model this is to introduce a time dependent lattice spacing a(t). 
This introduces a (pa/ a term in the equations of motion which, for a(t) increasing 
with time, leads to damping. Doing this, we found that each of the modes (now 
defined by spatially fourier transforming from comoving x to comoving k) red 
shifts with time. However, as long as the term does not dominate over the <fi 
term low momentum modes still become unstable and grow with respect to high 
momentum modes. In a real heavy ion collision, the expansion is anisotropic. 
The plasma expands much more rapidly along the beam direction than in the 
transverse directions [11]. While we have not attempted to model this explicitly, 
it probably means that correlated volumes do not grow large in the longitudinal 
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direction. The large expansion rate will rapidly damp modes of the pion field 
with wave vectors in the longitudinal direction. However, it seems possible that 
for a particular longitudinal position (or, equivalently [11], a particular rapidity) 
coherent fluctuations on length scales all the way up to the transverse extent of 
the plasma could develop. Another effect of the expansion is that after some time 
(perhaps of order 10 fm _1 ~ Im^ 1 [11]) the energy density is low enough that the 
description in terms of classical fields no longer makes sense. After this time, one 
has individual pions flying off towards the detector. From Figure la, it is plausible 
that at this time there are large long wavelength oscillations of the pion field. 

If the initial conditions or parameters are such that the initial value of the 
mean energy per unit volume E is much more than jv 4 no striking growth of low 
momentum modes occurs. Indeed when the total energy is large, the field will be 
able to climb far up the potential to \<f>\ 3> v, and therefore (0 2 )(£) will always be 
greater than v 2 , and no modes will become unstable. Requiring the energy not to 
be too large compared to the height of the potential at = can be seen either as 
a constraint on the parameters of the potential (V((j) = 0) ~ jv 4 ) for given initial 
conditions or on the initial conditions for a given potential. In the simulation of 
Figure 1, the total energy E ~ 2^t> 4 . As parameters (or initial conditions) are 
changed so that the total energy is increased, the growth of low momentum modes 
becomes less and less prominent, and when E ^ 10^f 4 , dramatic growth does not 
occur. 

As long as is centered around = initially, and as long as E is not too 
large, the choice of initial distributions has little qualitative effect. We have tried 
distributions ranging from = with chosen from a gaussian distribution to 
= with chosen from a uniform distribution with 2 < v 2 , and in all cases 
obtained results similar to those of Figure I. 

Let us attempt to estimate the conditions that might govern a realistic quench. 
We will make a crude estimate of the lattice spacing a and initial distributions for 
and using the one loop temperature dependent effective potential with H = 0. 
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At large temperatures, this is given by [15,16] 

y(0,T) = m 2 (T)0 2 + ^> 4 (10) 
where the temperature dependent mass is given by 

m 2 (T) = \{-v 2 + T 2 /4) . (11) 

While this is a crude approximation, it gives T c = 2v = 175 MeV which is coinci- 
dentally quite close to the lattice gauge theory value [17]. If we quench from an 
initial temperature Tj, the lattice spacing will be given by the correlation length 

at Tj 

a = m- 1 (T i ) (12) 

For Ti = 1.2 T c = 2.4 v, this gives a" 1 = 2.97 v = 260 MeV. 

We will model initial conditions with correlation length a by choosing <p (and 
(f>) independently at each lattice site. should therefore be distributed according 
to the probability distribution 

exp-^(V(T l )0 2 + ^ 4 ) . (13) 

Because it is cut off by the quartic term, this distribution has much less weight at 
large |0| than a gaussian. We approximate it by the truncated gaussian distribution 

exp-^m 2 (Tj)0 2 j for \<p\ < <p max ; for \<p\ > <p max (14) 

where (p ma x is chosen so that 2 has the same expectation value in the distributions 
(13) and (14). To choose initial conditions for <fi, we note that (0 2 ) = (0 2 ) in 
lattice units. (This is true for a quadratic potential, and is also true in general if 
equipartition applies.) Hence, we use the same distribution for as for 0. We now 
have a crude recipe for choosing initial conditions appropriate to model a quench 
from an initial temperature Tj to T = 0. 
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What does our crude recipe predict? From (13), we expect that in lattice units 

(fimax ~ 

(A/2) 1 /4 = 0.56. In fact, for T { = 1.2 T c , we obtain (pmax = 0.63 a = 
1.87 v. The value of <p ma x is crucial, because it determines over how much of 
the potential is distributed. We see that choosing gaussian distributed with 
variance v/2 as we did in Figure 1 was overly optimistic. For Ti = 1.2 T c , it is more 
reasonable to assume that is distributed according to (14) with (p ma x = 1-87 v. 
With these initial conditions, the energy per unit volume turns out to be E ~ 8^v 2 
and the pion field behaves as shown in Figure 3. Low momentum modes do indeed 
grow, but less dramatically than in Figure 1. Also, the system finds its way to a 
final state described by equipartition more quickly than in Figure 1. Nevertheless, 
there is a period of time after the quench when low momentum modes have much 
more power than higher momentum modes. 

If Ti is taken to be higher, 4> mSLX /v and 1/av increase, _E/(^t> 4 ) increases, and 
the phenomenon of interest is washed out. It is not unreasonable, however, to 
consider quenching from temperatures larger than but comparable to T c . If the 
initial temperature were much higher, the system could plausibly stay in thermal 
equilibrium until close to T c . Hence, choosing Ti = 1.2 T c and thus obtaining 
parameters and initial distributions like those of Figure 3 may be a reasonable if 
crude approximation to a real heavy ion collision. 

We have made many idealizations and approximations. Nevertheless, it seems 
possible to us that the essential qualitative feature of the phenomenon we have elu- 
cidated — long wavelength pion modes experiencing periods of negative mass 2 and 
consequent growth following a quench — could occur in real heavy ion collisions. 
Given the explicit symmetry breaking, one might have expected the dynamics fol- 
lowing a quench from a "generic" initial state to be featureless. The mechanism 
here discussed provides a robust counterexample that should be applicable in other 
contexts besides the QCD phase transition. Examples could include the reheating 
of some inflationary universe models and the quenching of spin systems which or- 
der at low temperatures, for which the interactions of the spins with phonons are 
less important than their interactions among themselves. If a heavy ion collision 
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is energetic enough that there is a central rapidity region of high energy density 
and low baryon number, and if such a region cools rapidly enough that it can be 
modelled as a quench, this will be detected by observing clusters of pions of sim- 
ilar rapidity in which the the fraction of neutral pions is fixed. This ratio will be 
different in different clusters and will follow a distribution like (2). We can think 
of no process besides a QCD phase transition with the chiral order parameter far 
out of equilibrium that could produce such a signature. 

Acknowledgements: We are grateful to David Spergel and Neil Turok for giving us 
the code created by the authors of [12, 13, and 14], and to both of them and Ue-Li 
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FIGURE CAPTIONS 



1) Time evolution of the power spectrum of the pion field and the sigma field 
at several spatial wavelengths after a quench at t — 0. and <p were chosen 
at time t — independently at each lattice site from gaussian distributions 
centered at the origin as described in the text. The parameters of the poten- 
tial were chosen so that for a lattice spacing a = (200 MeV)" 1 , the physical 
values TO-n- = 135 MeV, m a = 600 MeV, and = 92.5 MeV are obtained. 
The simulation was performed on a 64 3 lattice and the time step was a/10. 
After every two time steps, the spatial fourier transform of <f> was computed. 

— * 

For each component of 0, the power in all the modes with k = \k\ in bins 
of width 0.057a~ 1 were averaged. Figure la shows the time evolution of one 
component of the pion field. The curves plotted are (top to bottom in the 
figure) the average power in the modes in the momentum bins centered at 
ka = 0.20, 0.26, 0.31, 0.37, 0.48, 0.60, 0.71, 0.94, 1.16, 1.39, and 1.84. The 
initial power spectrum is white and all the curves start at t — at approx- 
imately 0.01 in lattice units. Hence, the longest wavelength pion modes are 
amplified by a factor of order 1000 relative to the shortest wavelength modes 
which are not amplified at all. Figure lb shows the time evolution of the 
sigma field. Only four modes are plotted — ka = 0.20, 0.48, 0.94, and 1.84. 
The vertical scale is different than in Figure la — the sigma modes do not 
grow nearly as much as the pion modes. 

2) Time evolution of the spatially averaged ((f) 2 ) for the same simulation whose 
results are shown in Figure 1. The horizontal line is at ((f) 2 ) = v 2 . 

3) Same as for Figure la, except that in this figure we have used a" 1 = 260 MeV 
and have chosen initial conditions from the distribution (14) following the 
crude recipe outlined in the text to model a quench from T{ — 1.2 T c to T — 0. 
Note that the vertical scale is different than in Figure la. All the curves 
start at about 0.03 at t — 0. With these more realistic initial conditions, 
long wavelength modes of the pion field grow less than in Figure la, but are 
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still significantly amplified relative to shorter wavelength modes. 
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